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ABSTRACT 

The management of endangered species as metapopulations is becoming increasingly common. Diverse as- 
pects of metapopulation dynamics and management have received attention in recent years. In particular, 
translocation of individuals between subpopulations of a metapopulation is practiced, or envisaged, for a 
variety of reasons and requires careful consideration. Linklater (2003) proposed that the number of individ- 
uals of each sex translocated into a target population for the purposes of maintaining genetic diversity could 
be chosen on the basis of parental investment theory. In this paper, following basic ideas in the parental 
investment literature, I propose a simple model to capture Linklater's proposal and provide a thorough 
mathematical analysis of the model. Granted the necessary species-specific biological information which 
would determine the model parameters in any instance of application, the analysis indicates that a practical 
algorithm can be constructed to generate the model's predictions for the optimal translocation. 
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1. INTRODUCTION 



Many endangered species now exist in fragmented habitat and are managed as meta-populations, the 
subspecies of black rhino Diceros bicornis of Africa providing notable examples. This situation poses many 
problems as regards reserve design, population dynamics, and management. In particular, reduced possibility, 
or even impossibility, of natural dispersal in such circumstances reduces genetic flow between subpopulations 
and may decrease productivity and increase risk of local extinction due to reduced fitness and stochastic 
effects (e.g., Soulc 1987). 

Translocation of individuals is one option for maintaining genetic flow between artificially isolated sub- 
populations of a given (sub)species. Given the inherent risks and practical challenges (e.g., Maguire et al. 
1987, Brett 1998, Miller et al. 1999, Fischer & Lindcnmayer 2000) involved in this practice, choosing strate- 
gies that optimize the reproductive return from translocation efforts and enhance metapopulation persistence 
is sound policy. A variety of issues have received theoretical treatment (the following list and literature are 
not intended to be exhaustive): which life stages arc optimal for translocation (Hearnc & Swart 1991, Robert 
et al. 2004); the role of captive populations (Tcnhmnbcrg et al. 2004, McPhce & Silverman 2004); decision 
making in the face of limited resources and uncertainties (Maguire 1986, Haight et al. 2000); metapopulation 
dynamics modelling (Hearne & Swart 1991, Wootton & Bell 1992, Lubow 1996); rate of genetic flow required 
(Wang 2004). 

As some of these studies emphasize, the goal of translocation is not just to merely supplement numbers 
but also to enhance reproductive performance and maintain genetic flow. Hence, the reproductive perfor- 
mance of the translocated individuals in the target population should be a primary concern and a measure 
of the success of a translocation. Linklater (2003) argued that managers planning a translocation can be 
viewed as parents, constrained by limited resources, investing in offspring and suggested that translocations 
be modelled in terms of parental investment theory. Application to a specific species will then depend on 
assumptions regarding reproductive behaviour and how individuals respond to investment. 

In a polygynous, sexually reproducing species, all females are expected to reproduce while only the 
fittest males are likely to obtain reproductive opportunities. Trivers and Willard (1973) argued that for 
such species natural selection should result in a mechanism that enables a mother to influence, presumably 
at conception, the sex of her offspring so as to gain the optimal return on her investment in that offspring 
from the resources available to her at that time. Speciflcally, they predicted that when a mother is in poor 
condition she should produce daughters rather than sons and invest more in the daughters she does produce, 
while when in good condition she should produce sons rather than daughters and invest more in the sons 
she produces. These predictions presume that the condition of the mother determines her ability to invest in 
her offspring, that investment in offspring translates into the condition of that offspring at sexual maturity, 
and reflect that, in the polygynous context of the considerations, males require a competitive edge to gain 
mating opportunities so only a mother in good condition is in a position to produce such a son while mothers 
in poor condition are better off investing their more limited resources in females. 

While there has been some controversy about what Trivers and Willard actually meant, see for example 
Carranza (2002), Cameron and Linklater (2002) have argued forcefully for the interpretation outlined above. 
Moreover, Cameron et al. (1999) and Cameron & Linklater (2000) found evidence for the Trivers- Willard 
(TW) effect in wild horses in New Zealand consistent with their interpretation. Although there have been 
many studies of the TW effect, with mixed results, Cameron (2004) conducted a meta-analysis of mammalian 
sex-ratio studies and reported that much of the inconsistency in results can be attributed to the variety of 
indices of maternal condition employed. Studies that utilized indices of maternal condition near conception 
provided almost unanimous support for the TW effect. In addition, Cameron (2004) described a possible 
mechanism for facultative adjustment of sex ratio, through the effect of blood-glucose levels on blastocysts, 
which would be consistent with the TW effect. Subtleties of life history of a species may confound the 
prediction of which sex beneflts more from maternal investment (Clark 1978, Leimar 1996, Sheldon & West 
2004) but does not necessarily undermine the general argument of Trivers & Willard. 

Linklater (2003) specifically proposed that the TW effect could be adapted as a model for how managers 
should invest resources in translocation operations for meta-population management of polygynous species 
such as black rhino; specifically, that the number of males and females to be released into a given target 
population, and the partition of resources amongst them, should be modelled as a function of the resources 
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available using the TW effect as the underlying prineiple. Optimization of sueh a model would then provide 
the recommended translocation strategy for managers to adopt in a given circumstance. The approach 
envisaged by Linklater would constitute an important component of metapopulation management. 

In this paper, I present a simple model of reproductive returns on investment in translocations. Lin- 
klater's specific proposal would be incorporated by appropriate choice of return functions for males and 
females. I describe the model in §§2-4, basic elements of which are adapted from ideas in Frank (1987). 
Solution of the model involves two optimization problems, which are solved in §§5 9. These solutions, along 
with numerical case studies, indicate both the generality and practicality of the model. 

2. INGREDIENTS OF THE MODEL 

The aim is to model the return on investment in a translocation, where, by return, I mean a measure of 
the reproductive success of the translocated animals (which reflects the transfer of alleles into the population). 

By fixed costs are meant necessary costs independent of the number of animals translocated which do 
not yield any return. In the modelling, the available resources shall mean the resources after the exclusion 
of fixed costs. 

Packaging costs are, by definition, costs per individual released which do not contribute to the animal's 
success, i.e., there is no return on packaging costs, but these costs must be included in the model. An 
example might be transportation costs, but in general packaging costs will be specific to each translocation 
effort and may well be gender specific. 

I assume that the resources invested in an individual yield a reproductive return that can be modelled 
by a suitable return function of said resources. I further assume that the return function is fixed for each 
gender in any given application. In reality, for a given species, a particular life stage may be optimal for 
translocation. One might therefore always choose individuals from that life stage or, with greater difficulty, 
allow different return functions for different life stages of each gender in the modelling. 

Let X denote the proportion of available resources to be invested in males and z the proportion invested 
in females, whence x + z = 1. I denote the return function for a male by m{x) and that for a female by 
f{z). The choice of return function for each gender models expected response to investment of resources and 
behavioural assumptions about the species; in particular, appropriate choices incorporate the TW effect into 
the model. Beyond the choice of return functions, there are three further steps in the modelling process: 
determining how to invest a given level of resources in a given gender; combining the results of the previous 
step to find the overall return from investment in both genders; computing the optimal overall return and 
hence the recommended translocation strategy. 

Note that even if there were an ideal return function that models the behaviour of individuals of a given 
gender of a given species in the context of unlimited resources, the return function employed in a given 
application must reflect the available resources. In particular, while returns no doubt level off in the face of 
unlimited resources, in the context of limited resources such behaviour may not be evident. This fact may 
mean that functions which are unrealistic return functions in the context of unlimited resources might serve 
over limited ranges of resources. It may also be the case that one knows only an approximation to the actual 
return function. 

For this reason, and for some mathematical generality, I treat a collection of candidate return functions 
which are nondecreasing and simply characterized by their rates of change, speciflcally, by their first and 
second derivatives. Each be regarded as approximations to the sigmoid in appropriate circumstances. 

3. INVESTING IN A SINGLE GENDER 

In this section I will employ the notation appropriate for males; it is modified for females in a straight- 
forward manner. When resources x are partitioned into x = Y^^=i with Xi applied to the i'th of n 
individuals, the total return is 

n 

M{x;n;xi,...,Xn) ■■='^m{x^). (3.1) 
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For given x, also define 

M{x) := max M{x;n;xi, . . . ,Xn)- (3.2) 

(n;2;i,...,a:„) 
xi-{ ^Xn—x 

A4{x) is the return from the optimal assignment of resources x to males with individual return function 
m{x). Similarly, one defines a function J^{z) for females. 

3. a Step-Function Returns 

Suppose each male requires a certain minimal investment a, beyond the packaging cost d, before any 
return on investment occurs and that any investment at and beyond a + d results in a fixed return r^. Then 
the return function for a single male is a step function: 

m(x) = (;■'"' + (3.a.l) 

^ ^ \0, if a; e [0,a + rf). ^ ^ 

By the nature of this individual return function, for resources x the optimal strategy is to release rim males 
where Um = [x/{a + d)]i, with [ ]j denoting the integral part of the bracketed number. It follows that 



M{x) = r„ 



a + d 



(3.a.2) 



3.b Linear Returns 

Suppose the returns are linear functions of investment beyond the packaging costs: 

™w={s'^-^' (3.b.i) 

If resources x are partitioned amongst n males, then 

n n 

M{x; n;xi,. . . , a;„) = m{xi) = a{xi — d) = a{x — nd). (3.b.2) 

i=l i=l 

Due to the nonzero packaging costs, the optimal strategy is to choose n = 1, whence 

M{x) = m{x). (3.6.3) 



3.C Increasing-Marginal-Return (IMR) Functions 

In this case the return function is increasing with increasing marginal returns (the graph is concave up), 
which may be characterized as follows: the function is zero on [0, rf] and on [d, 1] satisfies 

m{d) = m'{x) > m"{x) > 0. (3.C.1) 

One can use the result of the Appendix A to solve the optimization problem. As usual, it makes sense 
to assign each male invested in at least d, of the resources. Suppose one has so invested in n males. Put 
Uj = Xj — d > 0, j = 1, . . . ,n. Since m' is injective (m" > implies m' strictly increasing) on the relevant 
domain, (A. 5) has a unique solution yi = ■ ■ ■ = y^. But, as m" > 0, the Hessian in (A. 6) is positive 
definite at this critical point, i.e., one has a local minimum, not a maximum. So equal apportionment is the 
least optimal strategy in this case when resources are divided amongst several individuals. It follows from 
Appendix A that the optimal strategy is to invest in a single individual, whence M{x) = m{x). The linear 
return of 3.b is a limiting instance of this case. 
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3.d Diminishing-Marginal-Return (DMR) Functions 



In this case the return function is increasing with decreasing marginal returns (the graph is concave 
down), characterized by a return function m which is zero on [0,d], while on [d, 1], 

m{d) = m'{x) > m"{x) < 0. (S.d.l) 

One may require the function to have a horizontal asymptote, in which case limj^^oo m'{x) = 0. 

From Appendix A, since m" < on (rf, 1], m' is injective on this subinterval. Hence, for n satisfying 
(A. 2), (A. 5) again has a unique solution: j/i =••• = ?/„ is the only critical point and the Hessian (A. 6) is 
negative definite, i.e., this critical point is a local maximum. 

Thus, equal apportionment of resources x > d \s the optimal strategy amongst n individuals, for n 
satisfying (A. 2). For fixed x and n, the optimal return is therefore 

M{x;n) ■.= nm{x/n). (3.d.2) 

It remains to optimize with respect to n subject to (A. 2): 

M{x)= max M{x;n) (3.d.3) 

0<n<x/d 

with the understanding that when x < d there are no such n and A4{x) = 0, i.e., A4{x) = on [0,d]. Note 
that on {d, 2d\ , n = 1 is the only possibility, so in fact 

M{x) = m{x) on [0, 2d]. (3.d.4) 

Consider the function (j)x{v) := vm{x/v) for fixed x G {d, 1]. The restriction vd < x entails one may 
take the domain of (px to be {0,x/d). With v € {0,x/d), x/v € {d,oo), whence (j)x is smooth on its domain 
if m is smooth on (rf, oo), which is presumed. Differentiating with respect to v yields: 

(p'xiv) = m{x/v) - m'{x/v); 

(3.d.5) 

€iv) = ^m"ix/v). 

Since m" < on {d,oo), and x and v are positive, then < on its domain, i.e., 4>x is concave down on 
(0, x/d), whence has a local maximum which is its global maximum. Let Vx denote the value of v for which 
this maximum occurs, i.e., Vx solves (p'xi^x) = 0- Putting tx ■= x/vx G (d, oo), then 

'm'{tx) = txm'{tx); equivalently m^tx) = —■ (3.d.6) 

ix 

This equation has a simple geometric interpretation; namely, the tangent line to m at tx coincides with the 
line through the point {tx,Tn{tx)) and the origin (since they have a point in common and the same slope). 
As m is concave down for a; > d > and zero on [0, d], there is a unique line through the origin and tangent 
to m. This geometric fact characterizes tx and indicates it is independent of x. Thus, one may write t, 
defined by 

m'{t) = (3.d.7) 

If one did not require integral values, then the optimal n for given x would be = Vx = x/t, which is 
a simple linear function of x. Since, however, n must be integral, because of the concave down shape of (px, 
the optimal n for given x is the integer either side of x/t yielding the largest return: 



'X' 




'X' 




or 




.1. 


i 


.1. 



+ 1 according as 



> 



X 



/t]J \ltU J \[x/t] 



+ 1 TO 7 — ri 7 or not. 



(3.d.8) 
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Thus, for a given x there; are just two candidates for nx- Notice that for d < x < t, [x/t]i = 0, whence 
must be 1. In particular, ii t > 1, then Ux = 1 for all x G {d, 1]. 

3.d.9 Observation 

Given m as in (S.d.l), if the solution t of (3.d.7) satisfies t>l then 

M.{x) = m{x). 

li t < 1, however, as x increases, so does x/t and so the trend in Ux should be increasing. Note that 
although t is independent of x, (f)x is not; in particular, its critical point Vx = x/t changes with x, as does 
the domain {0,x/d). In general, A4{x) = nxm{x/nx)- 

Explicit examples of (S.d.l) are: on [d, 1]: 

h{x)=K{x-dy K>0, 0<s<l g{x) = K (^1 - (^^^ ^ s > 0. (3.d.l0) 

The functions h are derived from the cumulative distribution function for certain /?-distributions and were 
considered by Frank (1987). They increase without bound. The functions g exhibit a horizontal asymptote 
a.t y = K; consequently, g is bounded. Note also that for s » 1, g{l) « K, so by an appropriate choice of s 
one can ensure that g approaches its bound on the interval [0, 1]. 

For the function h, the solution of (3.d.7) is t = d/{l — s) > d but increases without bound as s — » 1. 
In particular, t>l s > 1 — d. For such s, (3.d.9) is valid. Note that as s — > 1, the function h is almost 
linear, whence the result just noted is consistent with the results for linear return functions. In general, 
Vx = {l- s){x/d). 

For the function g, the solution of (3.d.7) is i = \/s + ld > d. Taking the logarithm and using 
L'Hopital's rule shows that lims^oo \/s + 1 = 1 while limg^o y/s + 1 = e, i.e., d < t < de as s ranges from oo 
down towards 0. If \/r + ld = 1, for some r € (0, oo), then y/s + Id > 1 for < s < r. But t > 1 requires 
1 < de, i.e., d > 1/e w 0.37. In particular, if d < 1/e, then t < 1. So, if d is less than 0.37, then (3.d.9) is 
inapplicable. In general, Vx = [1/v^s + l\{x/d). 

3.e Sigmoid Return Functions 

One need not restrict to the logistic g{x) = 1/[1 + exp(— a;)]. One only requires, on [d, 1], 

m{d) =0 m'{x)>0 

(3.e.l) 

and for some p, d < p < 1, m"{x) 

Demanding a horizontal asymptote requires lima^^oo m'{x) = 0. Assuming that m is C'^, then 

m"{p) = 0. (3.e.2) 

For any given x, and n satisfying (A. 2), critical points of M{x;n,xi, . . . ,Xn) satisfy (A. 5). If xi, . . . ,Xn 
are all less than p, then as m' is injective on {d,p), as in (3.3) xi = ■ ■ ■ = Xn and such a critical point gives a 
local minimum of the return function. This phenomenon obviously occurs for x < p for which m is exactly 
as in (3.e). Assuming continuity of A^, for x < p, M{x) = m{x). 

When X > p, ii becomes possible to have solutions of (A. 5) with Xj ^ a:^, in which case Xj < p and 
Xk > p say. But the Hessian (A. 6) of such a critical point is indefinite. By (A. 5), m'{xj) = m'{xk)- To first 
order, diverting a small amount of the resources Xk from the fc'th individual to the j'th individual exactly 
balances. But, to second order, since m"{xk) < while m"{xj) > 0, the return from individual k will decline 
while that from individual j will increase. The overall change in return will depend, to second order, on the 
relative magnitudes of m"{xk) and m"{xj). Such critical points are saddle points. 



J > 0, for a; e [d,p); 
1 < 0, for X e (p, 1]. 
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If X is large enough that x/n > p. then x/n = xi = ■ ■ ■ = x,,, is a sohition of (A. 5) and (A. 6) indicates it 
yields a local maximum amongst all partitions of x into n amounts (including those partitions for which some 
Xj are less than p). Unlike in (3.d), a local maximum may fail here to be the absolute maximum because 
of the partitions containing Xj < p (the saddle-point critical points provide stationary points where the 
function may start to increase in certain directions and thereby ultimately exceed the value of the function 
at the local maximum). If so, the maximum for this value of n must occur on the boundary of the domain 
and the argument of Appendix A still entails that the optimal strategy for investing resources x is equal 
apportionment amongst n individuals for some n satisfying (A. 2). 

In summary, the optimal strategy when x > p is amongst the strategies nm{x/n), nd < x, x/n > p. 
Note that when x < p, there is no such n, whence n = 1 is the optimal strategy, so the earlier deduction 
for the case x < p may be included in this characterization. Hence, the situation may be regarded as a 
generalization of that in (3.d). 

Indeed, as in (3.d), one can consider the function (px{v) ■= vrn(x/v). The calculations (3.d.5) are 
unchanged. Again, one finds that there exists a t satisfying (3.d.7) which yields the unique critical point for 
(t)x- Put Vx '■— x/t as, before. In the present circumstances, 

€{v.) = Jm"(t). 

By the shape of the sigmoid graph, t > p, whence 'm"{t) < 0. Thus, Vx is the unique local maximum of (j)x- 
Moreover, since only those n for which x/n > p are of interest here, one can restrict v so that x/v > p, i.e., 
<^x is indeed concave down on the domain of interest. So, as in (3.d), one can solve (3.d.7) for t and then the 
optimal strategy is as in (3.d.8). In particular, x <t => nx — \ and (3.d.9) is valid for sigmoid returns too. 

Observe that \ip is close to 1, then the sigmoid function is approximately as in (3.3), while if p is close 
to 0, then the sigmoid is approximately as in (3.d). Moreover, the step-function return may be viewed as 
an extreme limiting form of the sigmoid, while the linear return is a limiting form of both (3.c) and (3.d) 
and thus an approximation to certain sigmoids. One would therefore expect some unifying perspective on 
the results obtained in this section. For example, the optimal strategy for step functions was similar to that 
in (3.d) being also equal apportionment but with the difference that it was amongst the maximum possible 
number of individuals. That said, the line through the origin that touches the nonzero part of the step 
function goes through the point x = a + d, i.e., for the step function, although (3.d.7) is not applicable, the 
geometric construction of t would give t = a + d, which is consistent with (3.d.7) in that t is the amount 
of resources one would ideally equally apportion (without the integral requirement imposed on v) amongst 
individuals. 

Indeed, a unifying viewpoint is as follows. For any return function m(x), the quantity m(x)/x is the 
return per unit resource. This quantity is maximal for the point of the graph of m(x) where the steepest 
line through the origin touches the graph. Denoting by t this value of x, m{t)/t is the maximal value of the 
return per unit resource. Since optimal invc;stment in the cases studied here never involves unequal partition 
of resources amongst several individuals, one may restrict attention to equal division of resources x into 
amounts y amongst individuals, which yields a return of {x/y).m{y) = x.ni{y)/y, and the optimal strategy 
would be ideally xm,{t)/t. Since the number of individuals must be integral, the actual optimal strategy is 
the best approximation to this ideal strategy. This argument reproduces, at least, the final result in each 
case (a)-(e) above. For step-function returns, t = a + d; for IMR functions (including linear), t = 1 and, 
although equal division of resources amongst several individuals is not actually optimal in any sense, the 
optimal strategy is to invest all resources a; in a single individual, which does amount to t = 1; for DMR 
and sigmoid functions, t is characterized by (3.d.7) (when m{x) is smooth, elementary calculus shows that 
m{x)/x has a critical point satisfying (3.d.7) which is a local maximum when m" < at the critical point). 

It is noteworthy that when one computes A4{x) = nxni{x/nx) for (3.d.l0) and the example sigmoid 
function of §6.c, the resulting graph oi M.{x) appears to represent a continuous function of x. I consider this 
issue in Appendix B, which results, for some return functions, in a simple procedure for obtaining M.. 

4. THE MODEL 
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In order for male and female reproductive returns to be comparable, male returns must be weighted 
by reproductive opportunities, which is achieved by multiplying by the appropriate sex ratio of the target 
population. Let S := (j)/ ji denote the ratio of reproductively active females to reproductively active males 
in the target population prior to the release and let Sp denote the same ratio after the release, where the 
subscript indicates that Sp depends on the actual partitioning of resources through the numbers of each 
gender released. While the interpretation of and may depend on the species considered, since the 
aim of any translocation is to introduce breeding stock into the target population, I will assume that each 
translocated animal is included in Sp. 

The model I pose is as follows. Let d denote the packaging costs for males and 5 that for females. With 
notation as in (3.1), define 

Ri{x) := max [M{x;n;xi, . . . ,Xn)Sp + F(l - x;k;zi, . . . , Zk)], (4.1) 

X\-\ \-Xn=X 

2lH |-Zfc = l — K 

where the maximization is over all permissible partitions: 

Ui := Xi — d G (0, X — nd); 
Wj := Zj — 5 & {0,1 — X — k6). 

Investment in a individual which yields no return is suboptimal for that gender. Moreover, such individuals 
should not count in Sp and so cannot enhance overall returns. Hence, such investment is suboptimal across 
genders and therefore need not be considered. 

For target populations with a large number of both males and females relative to the number of releases, 
5p Ri 5 and a simplified model results by replacing Sp by S. In this case, since S is constant, the optimization 
in (4.1) splits into separate optimizations of each summand, each of which then assumes the form as in (3.2) 
and (4.1) becomes 

R2{x):=M{x)S + J^{l-x). (4.2) 
One might wonder whether (4.1) ever reduces to 

R-i{x) :=M{x)S^+T{l-x), (4.3) 



n < 



k < 



1-x 



^ ^ .Tj — X, 



i=l 



i=i 



where Sx denotes the sex ratio of the post-release population for the optimal divisions of resources x amongst 
males and of 1 — a; amongst females. I shall consider this question case by case. I shall refer to (4.1) as 
Model 1, (4.2) as Model 2, and (4.3) as Model 3. 

It follows from the considerations in §3 that for each x € [0, 1] there are only finitely many possible 
strategies over which one must optimize to form M, T , i?2 and R\ , whence each of these functions is well 
defined. When m and / are continuous, it follows, as argued in Appendix B, that M. and T are each 
continuous, whence i?2 is continuous. Thus, Model 2 is well defined in that i?2 always possesses a maximal 
return on [0, 1]. The continuity of R\ depends, however, on the factor Sp, which may vary discontinuously. 
Thus, while bounded on [0,1], if discontinuous may not achieve its supremum (least upper bound). I 
will discuss this issue in the context in which it arises. 

In the case of (4.2), \i M.{x) and J^(l — .x) are smooth functions of x, then R!i{x) = Jv[' {x)S — T' {1 — x) 
and the critical points are given by 

M!{x)S = T'{l-x), (4.4) 

with a local maximum ensured by 

> R!i{x) = M"{x)S + r"{x). (4.5) 

The optimal strategy would then be at such a local maximum or at a; = or a; = 1. Since (4.4-5) require 
analytic expressions for M. and T, they are rarely a practical method of solution. 
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In the following sections. I c;ondiict case studies of Models 1 and 2. The purpose of these studies is to 
explore the mathematical behaviour of the model and present solutions to various scenarios. The biological 
relevance of these case studies is not explicitly considered, but see the discussion. For notational convenience, 
when (fe)males, are modelled with a return function as in 3a - e, I shall refer to step-function, linear. IMR, 
DMR, or sigmoid (fe)males, respectively. I do not explicitly study return functions of type 3c in the following 
case studies though it would not be difficult to do so. The results for linear returns indicate what to expect 
for returns of type 3c. 



5. STEP-FUNCTION MALES AND FEMALES 



The individual return function for males is as in (3.a.l), while the individual return function for females 



is: 



rf, ii X G [a + 6,1]; 
0, if a; e [0,a + 5). 



(5.1) 



Such return functions might be appropriate for modelling in situations where there is a regime of fixed 
investment per individual and only a knowledge of average returns modelled as a flat return; in particular, 
translocations called hard releases. 

As we saw in (3. a), the resulting optimal strategy for assigning resources x to males is 



M{x) 



a + d 



Similarly, the optimal strategy for assigning resources z to females is 



Hz) = rf 



(5.2) 



(5.3) 



Observe that ifa + d + a + 5 > 1, then the resources are sufficient only for investment in either males or 
females. If, further, both a + d < 1 and a + 6 < 1, then both options are available. For an investment 
exclusively in females the optimal return is ry.[l/(a + S)]i. For an investment exclusively in males, the 
possibilities are rmn4>/{n + n), n = 1, . . . , [l/(a + d)]j. Observing that n/(/z + n) is an increasing function 
of n, the optimal return is rm[i-/{(i + d)]iSp. (This last fact remains true if Sp is replaced by S.) Thus, one 
need only compare these two expressions to decide whether it is more optimal to invest in females or males. 

Now suppose {a + d) + {a + S) < 1, equivalently, a + d< 1 — a — 6, so that investment in both genders 
is possible. First consider Model 2: 



R2{x) 



a + d 



S + rf 



1 



a + S 



One can analyse this situation conceptually. Suppose 



1 = p{a + d) + ei, pa nonnegative integer, < ei < (a + d). 



(5.4) 



(5.5) 



On each subinterval of the form [s{a + d), (s + l)(a + d)), s — 0, . . . ,p—l, and on [p{a + d), l] , Al is constant 
and given by M{sia + d.)) while T is nonincreasing. Thus, R2 assumes its largest value on each of these 
subintervals at the initial point of the subinterval (though not necessarily uniquely so). It follows that the 
optimal strategy for R2 occurs at one of: 



s{a + d), 



s = 0, 



(5.6) 



Note that there may be other values of x at which an equally optimal strategy occurs. Such values arise 
because neither A4 nor may change for small changes in x. For example, the analogous approach to 
(5.7) based on increments of a + S would proceed in such increments down from 1 towards 0. Note that 
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the resulting points y = 1 — (.{a + (5) may not have any points in common with (5.7). If not. this fact is 
explained by the previous remark. As long as one finds some x that gives the optimal strategy, it is a matter 
of convenience what form one takes; unless there are restrictions on the number of males or females available 
for translocation. Here I take the form given by (5.7); but the second approach will prove convenient in the 
next section for case studies with step-function females and male returns of other forms. 
If a + > a + (5, write 



a-V d = q{a + d) + 62, q & nonnegative integer, < 62 < a + 5. 
To compute the difference in return R2 between x = s{a + d) and x = {s + l){a + d), suppose 
1 — s{a + d) = is{a + 6) + Cs, is a, nonnegative integer, < < a + 5. 

Then 



(5.7) 



(5.8) 



while 



R2{s{a + d)) = STmS + rfis 

'1- {s + l)(a + d) 



R2{{s + l){a + d)) = {s + l)rmS + rf 



S 



By assumption < 1 — (s + l)(a + d), whence 1 — {s + l){a + d) + {a + d) = 1 — s{a + d) yields is+i{a + 6) + 
Cs+i + q{a + 5) + €2 = £s{a + 6) + es, which entails £s > 9- Hence, 1 — {s + l){a + d) = {is — q){a + S) + es — 62, 
and 



1- {s + l)ia + d) 



Thus, 



I.e., 



a + S 

R2{{s + l){a + d)) - R2{s{a + d)) =rmS 



£s-q, ifes>e2; 
£s-q-l, if < €2- 



R2{{s + l){a + d))>R2{s{a + d)) <^ 



Hence, with a + d > a + S, when 



> 



a + d 
a + S 



> < 



1, 



r/g, if es > €2: 

rf{q+ 1), if < £2, 



if eg > £2; 
1, if < £2, 



















a+<5 



(5.9a) 



it does not diminish returns to divert a further a + d resources from females to males, i.e., (5.9a) implies the 
optimal return occurs at a; = p{a + d), i.e., invest in as many males as possible and the residual 1 — p{a + d) 
in females if possible. If, on the other hand. 



< 



a + d 



a - 



(5.96) 



it never pays to divert a+rf from females to males, so the optimal return is obtained by investing all resources 
in females. Any left over resources will be less than a + S, whence less than a + d, and therefore insufficient 
for any male. Finally, if 



a - 



1>^> 



a + d 



a - 



(5.9c) 



then at any given x = s{a + d) it may or may not pay to divert a further a + din resources from females to 
males. In this case, one must examine each possibility in (5.6). 

Now suppose instead that a + d < a + S. One can analyse this case by suitably interchanging the roles 
of males and females. Let A4' be the function (5.2) that results from the female model parameters a, 6, rf 
and J^' the function (5.3) that results from the male model parameters a, d, rm- Then A4'{x) = T{x) and 
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T' (z) = M.[z). Consider j\A'{y)(\l S) ^ T'(\ — y). Since is a constant, optimization of this expression wrt 
y is equivalent to optimizing M.'{y) + — y)S wrt y, i.e., equivalent to optimizing T{y) + M.{1 — y)S 
wrt y, which is equivalent to optimizing M{x)S + !F{1 — x) wrt x = 1 — y. Hence, interchanging the roles 
of males and females and inverting the sex ratio yields a problem that is equivalent to the original model. 
Thus, one solves the transformed model as in the previous paragraph; if the solution is to apply resources y 
to 'males' represented by M' , then the actual optimal solution is to apply resources y to actual females. 

For a direct approach to the case when a + d < a + 5, the question is whether it pays, at a; = s{a + d), 
to divert a further a + d in resources from females to another male. The gain from the additional male is 
TmS. Again write 1 — s{a + d) = ^s(a + (5) + Cg. Since a + d < a + 5, there is no loss in return from females if 
a -\- d < eg and a loss of if a + > (one loses at most one female in this case to gain a + for a male) . 
In the latter case, the net change is rmS — rf. Whether the latter case occurs depends on e^, which depends 
on the model parameters. In any case, one can assert that if 

a + d<a + S and — >1 (5.10a) 

r/ 

then it does not diminish returns i?2 to divert further resources from females to males, whence the optimal 
return is obtained by investing in as many males as possible (the residue is too small to invest in a female). 
On the other hand, if 

r S 

a + d<a + d and < 1 (5.106) 

rf 

then it will not pay to divert a + d to another male if doing so reduces the number of females invested in 
(by one). In this case, one can only say the optimal return occurs at some point of (5.6); computation will 
determine which. 

Observe that if one solves this scenario using the symmetry argument, then (5.9b) asserts that if 
Tf/vmS < [{a + 5) / (a + d)]i then the optimal strategy is to invest in all "females", i.e., in actual males. Since 
a + d < a + S, then [(a + S)/{a + d)]i > 1. Thus, (5.10a) is consistent with (5.9b). If (5.10b) holds, any of 
(5.9a-c) may in principle hold. 

Note that the conditions (5.9 10) only involve the ratio of VmS to rf and the ratio of a + d to a + S. 

It is a simple matter to compute i?2 for specified values of x, in particular the values (5.6), with the 
model parameters as variables (for example in MS Excel) so that one can easily numerically solve the model 
and generate examples confirming (5.9-10). For example, with a — 0.05, S = 0.1, a = 0.2, d = 0.05, rf — 1, 
rm = 1-5, S = 1, then 1 < (a + d){a + 6) and [{a + d){a + S)]i < rm/rf < [{a + d){a + S)]i + 1. Consistent 
with (5.9c), one finds the optimal strategy occurs at x = 0.25, i.e., s = 1 in (5.6), so a mixed strategy of 
investing in 1 male and 5 females yields the optimal return for the given model parameters. 



Now consider Model 1 under the assumption {a + d) + {a + S) < 1; equivalently, a + d < 1 
Then, Ri{x) is obtained by maximizing 



rrr,n- 



fj. + n 



rfk, 



with respect to n and k subject to 



< n < 



a + d 



0<k< 



a + 5 



5. 



(5.11a) 



(5.116) 



but where n = only on [0, a + d) and k = only on (1 — (a + S), 1]. As a function of n, the first summand 
is increasing and so is maximal for n = [x/{a + d)]i. With this choice, (5.11a) is also an increasing function 
of k, whence maximal for fc = [(1 — x)/{a + S)]i. Hence, 



Ri{x) = M{x)S^ + J^{1 - x) = R3{x), 



(5.12) 



i.e., (4.1) reduces to (4.3); Model 1 reduces to Model 3. 

So consider Model 3. Recalling (5.5), on a subinterval [s{a + d),{s + l){a + d)), s = 0,...,p— 1 or 
[p{a + d),l\, since 

-I I III ~t- 

X 



Raix) = r„ 



a + d 




1 



a + 5 
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then as [(1 — x)/{a + 5)\i is nonincrcasing and [x/{a + d)]i is constant, then Ad{x) is c;onstant, T{1 — x) is 
nonincreasing, and Sx is also nonincreasing (the denominator is constant and the numerator nonincreasing) . 
Hence, Rz{x) takes its maximum value on each of these subintervals at the initial point of the subinterval 
(though not necessarily uniquely). Hence, as for _R2, one can find an optimal strategy from amongst the 
possibilities (5.6). This observation limits the numerical analysis required. Note, however, that the gain in 
diverting a + d from females to another male is offset not only by the loss of the return from the females but 
also by the decrease in Sx through the numerator, so that (5.9) is not valid for Model 1. 

Suppose a = 0.05, 5 = 0.05, a = 0.2, d = 0.025, r™ = 4, r/ = 1, ^ = 10, /x = 5 (whence S = 2). Then 
(a + d)/{a + 6) > 1, Srm/rf = 8 > [(a + d)/{a + S)]i + 1 = 3, so in Model 2 (5.9a) predicts the optimal 
strategy would be to invest in as many males as possible. Indeed, for Model 2, x = 0.9, i.e., s = 4 in (5.6), 
is the optimal strategy, investing in 4 males and 1 female for a return of 33. For Model 1, however, one 
computes that the optimal return is 22.5 and can be achieved with s = 3 in (5.6) to give x = 0.675, i.e., by 
investing in 3 males and 3 females. In Model 1, the strategy x = 0.9 yields a return of 20.5. 

6. STEP-FUNCTION FEMALES AND NON-STEP-FUNCTIONI 

MALES 

In this section I assume the individual female return function is as in (5.1) while the individual male 
return function is assumed to be a nondecreasing function (other than a step function) with a packaging 
cost of d. I restrict attention to the case for which d+ {a + 5) < 1, equivalently d < 1 — {a + S), so that 
investment in both genders is possible. For Model 2, one therefore has (4.2) with J^{z) as in (5.3). Write 

1 = n{a + S) + e, n a nonnegative integer, < e < a + 5 (6.1) 

and choose the nonnegative integer k so that 

w:=l-k{a + 6)> d but 1 - {k + l)ia + S) < d, (6.2) 

i.e., subtracting increments of a + S from 1, w is the last quantity so obtained that is greater than d. Observe 
that, on each of the subintervals 

(1 - (a + S), 1] , (1 - 2(q + 5),l-{a + 5)] , {d, l-k{a + 5)] , (6.3) 

the function ^(1 — x) is constant, given by its value at the right-hand endpoint. 

6.4 Lemma 

Whenever m{x) is a nondecreasing function of x, so is A4{x). 

Proof: For resources y > x, one can of course obtain the return A4{x) with investment of resources 
X, leaving y — x further to invest, so one can do no worse than ^A{x), i.e., ^A{y) > J^{x). For m as in 
(3.d) or (3.e), for example, this result is explicitly of the form: M{y) = nym{y/ny), which is at least as 
large as nxm(y/nx) since Uy is by definition the optimal equitable division of resources y; but since m is 
nondecreasing, nxm{y/nx) > nx'm{x/nx) = M.{x). 

It follows from (6.4) and the constancy of !F{1 — x) on each of the subintervals of (6.3) that R2{x) is 
nondecreasing on each of these subintervals, whence its maximal value on each subinterval of (6.3) occurs 
at the right-hand endpoint. Moreover, since R-iix) = F{1 — x) on [0,(i], then R2{x) is nonincreasing on 
[0,d]. Thus, the optimal return in Model 2, whenever females have step-function return, occurs at one of the 
points: 

1, I -{a + 5), l-2{a + 8),...,l- k{a + 5)=w, 0. (6.5) 

For Model 2, it therefore suffices to compute the return R2 for this finite list of possibilities and select the 
maximum value to determine the optimal return. 
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Now consider Model 1. From (4.1) and (5.1), it is clear that the optimal strategy requires investment 
of 1 — a; in := [(1 — x)/{a + 5)]^ females so that Ri takes the form 

Ri{x) = P{x, nx){(f> + ki-x) + Tfki-x, (6.6) 

where P{x,n) := nm{x / n) / {fi + n) and Ux is the nonnegative integer maximizing P{x,n) subject to (A. 2) 
{ux = 1 for 3.b-c). Hence, if y > a;, because P{y,ny) > P{y,n) for any n satisfying (A. 2), 

_ nym{y/ny) nxm{y/nx) nxm{x/nx) _ , , 
ji + Uy IJi + Ux At + rix 

Thus, P{x,nx) is a nondecreasing function of x and the remaining terms of (6.6) arc constant on the 
subintervals (6.3), so R\{x) is nondecreasing on these subintervals and thus takes its maximal value on each 
subinterval at the right-hand endpoint. As R\{x) = J^{1 — x) on [0, d], then as with Model 2, it suffices to 
compute the return i?i for the finite list of possibilities (6.5) and select the maximum value to determine 
the optimal return. 

I present details of specific cases in the following subsections. 

6. a Linear Males 

The individual return function for males is given by (S.b.l). As usual, I suppose d + {a + 6) < 1, 
equivalently, d < 1 — a — 6, in order that investment in both genders is possible. 
For Model 2, one has 



R2ix) = < 



a+5 



if X e [0,d\; 

(6.a.l) 

a{x-d)S + rf , ifxG[d,l]. 



I will derive a complete characterization of the optimal strategy in terms of the model parameters. 
Assume that 

a{a + S)S>rf, (6.a.2) 

i.e., the marginal return on investing a+d in a male (beyond the initial investment of d) exceeds the marginal 
return from investment of the same resources in a single female. Then, 

i?2(l) > i?2(l - a - (5) > • • • > i?2(l - (A: - l)(a + S)) > i?2(w) (6.a.3) 

since each term is the previous term plus —a(a + S)S + rj. Thus, granted (6. a. 2), (6.5) implies that i?2 has 

its maximum value on (d, 1] at a; = 1. It then follows from (6.5) that the larger of -^2(1) and i?2(0) is the 
maximum value of i?2- 
If, on the other hand, 

a{a + 5)Sp < r/, (6.a.4) 

then 

i?2(l) < i?2(l -a-5) <■■■ <R2{l-{k- l){a + 8)) < R2{w), (6.a.5) 

whence from (6.5) the maximum value of R2 on {d, 1] occurs a.t x = w. Hence, by (6.5), the maximum value 

of i?2 occurs at X = w or X = (the return on investing a + S in another female is greater than the marginal 
return on the same investment in a male, so the optimal strategy is to invest in as many females as possible, 
either exclusively females or the maximum number with enough left over to invest in a male with return). 
Clearly, if i?2(0) < -R2(l), then as -R2(l) < R2{w), then the optimal strategy occurs at x = w. 

To determine the optimal strategy, otherwise, i.e., when (6. a. 4) holds and i?2(0) > i?2(l), recall (6.1), 
and suppose first that d < e, whence the k of (6.2) equals n and w = e. In fact, irrespective of the relative 
magnitudes of -R2(0) and .^2(1), when d < e, 

R2{0) = TfU < rfn + a{e - d)S = i?2(e) = R2{w), 



13 



i.e., X = is not optimal. 

Finally, suppose instead that d> e, whence the k in (6.2) is less than n. Now, using (6.1), 



R2{w) = rfk + a(l - k{a + S)~d)S 

= rfk + a(l - d)S - ak{a + S)S 

< rfk + a(l - e)S - ak{a + 5)S 

= rfk + an{a + 5)S - ak{a + 5)S by (6.1) 

= rfk + a{n - k){a + 6)S 

<rfk+{n~k)rf by (6. a. 4) 

= r/n = E2(0). 



In summary, the optimal return must occur at a; = 1, a; = 0, or the interior solution w given by (6.2), 
and these possibilities occur as follows: 

i?2(l) > -R2(0) & aS{a + 6)>rf then x = 1 
(if ii2(l) = -^2(0), then a; = and a; = 1 give the same return); 
R2{0)> R2{1) kaS{a + 6)>rf then a; = 0; 
R2{1) > R2{0) k aS{a + 5) <rf then x = w; 



(6.a.6) 



Each of these scenarios is realizable with appropriate choices of model parameters. 

Now consider Model 1. Investing in more than a single male decreases the return from males both 
through packaging costs (as in 3.b) and through the weighting factor Sp. Hence, Model 1 reduces to Model 
3 and (6.6) becomes 



Ri{x) = < 



a+S 



l-x 



a+S 



if a: e [0,c/]; 
if a; e [d, 1]; 



(6.a.7) 



which is in fact continuous at a; = d. 

The optimal return must occur at one of (6.5) but the varying value of Sp now complicates matters. 
One has: 



i?i(l) =a{l-d) 



M + 1 



Ri{l- {a + 5)) =rf + a{l- {a + 8)- d)'^—- 



a(l - d) 



+ 



r/ — a{a + 5) 



Ai + 1 



i?i(l-2(a + (5)) =2r/ + a(l-2(a + (5)-d)^^ = r/ + a(l - (a + 5) - d)^^ 



+ 



r/ — a{a + 5) 



/U + l_ 

(6.a.7) 



Ri{w) =i?i(l-fc(a + (5)) = krf + a{l - k{a + 5) - d) 



<t> + k 
At+ 1 



[k - l)r/ + a(l - (fc - l)(a + 5) - d) 



M+ 1 
i?i(0) = r/n. 



+ 



r/ — a{a + 5) 



/U + 1 
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For r = 1, . . . , k, consider the final expression for i?i (l — r{a + S)) . The first bracketed quantity is larger than 
i?i(l — (r — l)(a + S)). The second summand can be either positive or negative. Thus, depending on these 
values, it seems likely that the optimal return could occur at any one of (6.5), yielding a more complicated 
situation than for i?2- 

Consider two possible simplifications. Suppose 

Tf > a(a + 6)- > a(a + 6)- ^— — - > ■ ■ ■ > a(a + 6)- — -, (6.a.8) 

noting that the only assumption here is the first inequality. Then from (6. a. 7) one sees that 

Ri{w) > Ri{l - {k - l){a + 5)) > ■ ■ ■ > Ri{l - {a + S)) > Ri{l), (6.a.9) 

whence the maximum of i?i on (d, 1] occurs at x = w. It follows that the maximum of Ri occurs at x = w or 
a; = if (6. a. 8) holds. This result is an analogue of the results obtained above for R2 under the assumption 
that r/ > a{a + d)S. 

If on the other hand 

rf<a{a + 5)^<a{a + 6)^<---<a{a + 6)^, (6.a.9) 

where again only the first inequality is an assumption, the second summand in each line of (6. a. 7) is now 
negative, but whether sufficiently negative to make i?i(l — r{a + 5)) smaller than Ri(l — (r — l)(a + (5)) 
is undetermined, which reinforces the possibility that any of (6.5) might give the optimal return, depending 
on the values the model parameters take. 

In any event, one need only examine the possibilities in (6.5) to find the optimal strategy. 

Consider a = 0.05, 5 = 0.05, r/ = 1.5, a = 10, d = 0.05, 5 = 3 whence (6.a.6) is satisfied. With 
Sp = 30/10 there is no difference between Model 2 and 1 as regards optimal strategy: choose x = 1, i.e., 
invest all resources in a single male. But, in Model 1, putting Sp = 6/2, which involves no change for Model 
2, the optimal strategy for Model 1 is now: choose a; = 0.6 = 1 — 4(q: + d), which gives a return of 24.3, by 
investing only 60% of resources in the male and the remaining 40% equitably amongst 4 females, versus a 
return of 19 for x = 1. Moreover, x = 0.6 is not a possible optimal strategy for Model 2 at all {w = 0.1). 
This example underscores the difference between Models 1 and 2. 

6.b DMR Males 

With / as in (5.1) and m as in (3.d.l), one can write for Model 2 

(6.b.l) 

where, as usual, is the optimal number of males in which to invest resources x, with 71^ = understood 
for X G [0,d\ so the first term reduces to zero on [Q,d]. Assume that d + a + 5 < 1 so that investment in both 
genders is possible. 

Suppose TO is known explicitly. If one can solve (3.d.7) cxplicity for t, n,,. can be determined from (3.(1.8): 
otherwise, one can approximate t with sufHcient accuracy for numerical computations. For certain return 
functions, one can compute M from the observations in Appendix B. Alternatively, one can compute the 
quantity nm(x/n) for the allowed values (A. 2) of n and determine rix by inspection. Thus, for an explicit 
m, one can compute M, and then i?2, at each of the values in (6.5) and numerically determine the optimal 
return in Model 2. Can one say more? 

If one employs the approximation 

Tlx = |, (6.b.2) 



-^2(2;) = riajTO ( — 1 S* ■ 



a + S 
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with t the sohition of (3.d.7), then in effect the male return function becomes piecewise-linear and one obtains 
an analytic expression for the resulting approximation to R2: 



{xm{t)/t)S + rf , iix€{d,l]. 

This expression is of the form M'^{x)S + !F{1 — x) with 

0, if X e [0, d 



A^"(a;) := "',.,, " T (6.b.4) 

^ ' yxm{t)/t = xm yt), itxS(a, IJ; ^ ' 

noting (3.d.7). Thus, (6.b.3) differs only from the form of (6.a.l) in that M.°'{x) ex x rather than {x — d), 
which makes A^" discontinuous at a; = (note that the packaging costs d have the effect of making the 
approximation for male returns piecewise linear). Working through the argument of §6. a as applied to 
(6.b.3), one obtains a modified version of (6. a. 6): 

^2(1) > ^2(0) & m'{t}{a + 6)S > rj then x = 1 
(if i?f (1) = -R|(0), then x = and x = 1 give the same return); 

R^{l)>R^{0)&Lm'(t)ia + 6)S<rf then x = w; (6.b.5) 



i?^(0) > R^{1) & m'{t){a + S)S <rf k. 



d < e, then x = w; 
d> e, then x = or 



Because i?|(l) = m'{t)S while R^{0) = rf [l/{a + 6)]^ < rf/{a+6), rf < m'{t){a+d)S entails R^{0) < i?f (1) 
and the second possibility in (6.a.6) will not in fact occur here. But, while 

a + a + 

m'{t)S > i?2(l) so in the actual model i?2(l) < -^2(0) may not be ruled out when m'{t){a + 5)S > rf. 
Moreover, in place of (6. a. 4), one only finds R^{w) < -^2(0) + rn'{t)Se, so the final case in (6. a. 6), is here 
left unresolved, the condition d>e being uninformative; one must simply compute the two possibilities. 

Note that this approximation, if useful at all, is useful for obtaining the value of x for which i?i or R2 is 
maximized. The optimal strategy must still be interpreted as equitable partition of resources x in males, 
with Tlx given by (3.d.8) (and equitable partition of 1 — x in females of course). If the optimal return occurs 
in this approximation at x = 1, then of course Ux = 

How useful is the approximation (6.b.2) and its consequence (6.b.5)? The error in (6.b.2) is less than 
one, i.e., \nx — x/t\ < 1. Since Vx '■= x/t is the actual value at which Xx{v) — vm{x/v) achieves its maximum 
(§3.d), then the approximation (6.b.2) overestimates the actual return on the investment in males. The error 
will therefore be magnified by S > 1 and diminished by S* < 1. Thus, when (6.b.5) favours investment in 
females over males, that deduction should be robust. When (6.b.5) favours investment in males over females, 
that deduction may be questionable if the overestimation of male returns is significant. 

In various trials with m = 5 as in (3.d.l0) (especially with s = 1 for which Xx{v) is quadratic) the linear 
approximation (6.b.2-5) successfully predicted the optimal strategy. But, with d = 0.07, K = 5 and s = 1 
in (3.d.l0) and with S = 0.04, a = 0.03, and r/ = 1.4 in (5.1), then m'(t)(a + S)S < rj and R^{1) < R^{0) 
so (6.b.5) implies the optimal return occurs a.t x = w = 0.02 or x = 0. In fact, i?2 has its maximum 
value of 19.81 at x = li; (though R2{0) = 19.6, only sHghtly less) while R2{w) = 19.31 which is less than 
-^2(0) = 19.6, the actual maximal return being i?2(0.16) = 19.61. So, in this case, the approximation is 
slightly misleading. Note that 0.16 = 1 — 12(a + S), so this result is consistent with (6.5). It may well be 
that the difference in returns between the strategy suggested by (6.b.5) and the true optimal strategy, as in 
the example just presented, is of no practical significance. For an explicit choice of function for m, one can 
compute the error involved in utilizing R^, but since the optimal return for R2 can be obtained readily from 
(6.5) the approximation (6.b.2-5) is mostly of theoretical interest in suggesting that Model 2 with DMR 
males behaves similarly to Model 2 with linear males. 
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Now consider Model 1. Since for a fixed number of males invested in, the optimal partition of resources 
is equitable distribution, for any x £ [0, 1], the expression to be optimized to obtain Ri is 



nm{x/n) 



+ rfk 



< n < 



< fc < 



a - 



(6.b.6) 



where n = on, and only on, [0, rf] and fc = on, and only on, {1 — a — 5,1]. It is clear that k should take 
its maximal value. On (d, 1], therefore, (6.b.6) becomes 



nm{x/n) 





a+S 


:\ 




'l-x' 


\ n + n 





< n < ^, 
a 



(6.b.7) 



in agreement with (6.6), and where for x € [0, d] it is understood that the first summand vanishes. 

The final optimization required to obtain Ri is with respect to n. Analogous to the procedure employed 
to obtain in §3.4, for x e {d, 1], consider the function 



^x{v) ■= vm{x/v) 



(6.b.8) 



defined on the interval (0, x/d), where ~ (f) + [{1 — x)/{a + 5)]i. We aim to find the maximal value of ip^ 
as a function of v for each x G {d, 1] and deduce therefrom the optimal n^,. By direct calculation, and after 
some simplification: 



m 



III/ \ '^kx ii / \ , II ( ^\ 

which is not as simple as was obtained for (p^ in §3.4. First note that tp^ is zero when 

\V J jJL + V \V J \V J 

If Vx is a solution of this equation, with Tx ■= x/vx, then (6.b.ll) becomes 



m{Tx) 



M + ^/tx 



m'{Tx). 



(6.b.9) 
(a.b.lO) 

(e.b.ll) 
(6.6.12) 



Interpreting this equation geometrically, since 12/(11 + x/tx) < 1, is a point where the line through the 
origin and the point {Tx,m{Tx)) is steeper than the tangent line to m through the same point. Prom the 
shape of m, it is then clear that Tx > t. From (6.b.l0), %p'x{vx) < since m"{Tx) < 0, i.e., each stationary 
point of %px is indeed a relative maximum. Observe, from (6.b.l0), that 'ipx{v) is negative wherever 'ip'xi'^) > 0) 
but that when i)'x{v) < 0, Vx(^) ™=^y become positive. In particular, assuming V'x is C^, as v increases away 
from Vx one expects ipl^ may eventually change sign. But since there are no relative minima, a change in 
concavity cannot result in ^x exceeding ^px{vx) whence there can only be the one relative maximum. At 
V = x/d observe that (6.b.8) vanishes, whence, using (6. b. 8-10), i/'a can be extended to {Q,x/d\ in a 
fashion. Since one is only interested in the restriction of ipx to [1, x/d], on which tjjx is boimdcd and positive, 
and ipx has a unique relative maximum on {Q,x/d) and no other critical points, its maximum on [l,a;/(i] 
must be at Vx, whence Ux must be one of the integers either side of (the unique solution to (6.b.ll)) Vx- 
Thus, one docs obtain a result analogous to (3.d.8). 

But, we do not have an explicit characterization of Tx, only the implicit characterization of (6.b.l2) and 
we do not know the dependence of Tx upon x. To explore these issues, rewrite (6.b.l2) in the form 



m{Tx) - m'{Tx)Tx _ x_ 
m'{Tx) IJ,' 



(6.b.l3) 
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Wli«i /i S> 1. the right-hand side is approximately zero, whenee the numerator of the left-hand side is 
approximately zero, i.e., Tx ~ t. Note that x/ n is a simple linear increasing function of x. The expression 
(to(z) — m'{z)z)/m'{z), as a function of z, has derivative —mm"/{m')^, which is indeed positive. Thus, as 
X increases, the left-hand side of (6.b.l3) must increase, and since its a strictly increasing function of Tx, it 
must be the case that Tx increases with x. Since Tx > t, then as x increases, m'{Tx) decreases more rapidly 
than m{Tx)/Tx- It follows from (6.b.l2), that Vx = x/tx increases with x. 
In summary, one finds 



Ri{x) = nxm{x/nx) 





l-x 
a+S 


.] 


+ rf 


'l-x' 


\ l^ + rix 


a + S_ 



(6.b.l4) 



where, for x € (rf, 1], Ux is the integer either side of Vx — x/tx that maximizes the first summand in (6.b.l3), 
while Ux = iov X £ [0,d\. We have seen that Tx > t and that both Tx and Vx are increasing functions of 
X. Since Tx > t, then Vx ■= x/tx < x/t, i.e., not surprisingly, in Model 1, Vx, whence possibly n^, is smaller 
than in Model 2, since the value of males is discounted by their own numbers. Note that Model 1 has not 
reduced to Model 3. In any event, one computes R\ for each of the possibilities in (6.5) to determine the 
optimal return; for each such x, if necessary, one can compute nm{x/n)/{p. + n) for the allowed values (A. 2) 
of n and thereby determine Ux , thence Ri [x) . 

Turning to the explicit examples of (3.d.l0), for h one finds that (6.b.l3) becomes [{1 — s)Tx — d]/ s = x/ fi, 

i.e., 



sx + dfj, 



= t- 



sx 



As expected, Tx is an increasing function of x and for ^ 1, r^; ~ Furthermore, 

— s) 



X 
Tx 



(6.b.l5) 



(6.b.l6) 



which is an increasing function of x with horizontal asymptote at y = — s)/s. But a; = 1 is the upper 
limit for our interests, hence Vx < /^(l — s)/(s -|- d/u). Thus, for h, one has a fairly simple situation. 
Now consider g. For (6.bl3), one obtains 



t|+1 - (s + l)d'Tx ^ X 

sd^ jj. 
which can be rewritten as the polynomial equation 

p{z) := z'+^ -{s + l)d''z - {x/ii)sd' 

Observe that 



0. 



(6.b.l7) 



(6.b.l8) 



p{t) = p{</I+ld) = (s + l)</7+ld'+^ - (s + l)</7TTd'+^ - {x/iJi)sd' = — sd' < 0. 

A* 

By Descartes's rule of signs, (6.b.l8) has at most one real positive root. Since p is positive for large z and 
negative a.t z = t, then p has exactly one positive real root, somewhere on (t,oo). Also note. 



p{l) = l-{s + l)d' - -sd' <1. 

A* 

So, if p(l) > 0, then t < Tx < 1- Otherwise, Tx > 1, in which case, as in (3.d.9), Ux = 1- 

Now T (I drop the subscript for convenience) satisfies (6.b.l8). Differentiating this equation implicitly 
with respect to x, one finds 

s + 1 /i — a* 
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and as r > (i then > so r' is positive, confirming r is an increasing function of x. As Vx = x/t, 
differentiating with respect to x yields, 

< = (6-b-20) 
Substituting (6.b.l9) into (6.b.20), the numerator is 

, (s + 1)t«+i - (s + l)d^T - sxd^/n sr^'+i 

J- — XT -— — - - - ■ = 

(s + l)(r^-rf«) (s + l)(r*-rf«) 

since r satisfies (6.b.l8). As — rf* > 0, it follows that 

[Hi 
{s + 1){t'' - d") 



< = , fw_. > (6.b.21) 



confirming that Vx is an increasing function of x. In general, it does not appear possible to give an explicit 
formula for the positive real root of (6.b.l8), i.e., a formula for r as a function of x. 
Consider the particular case s = 1, however. Then (6.b.l8) is the quadratic 

— 2dz — xd/ n = 



with solutions d ± y/d"^ + xd/ ji. The choice of minus sign clearly gives a negative number so 

Tx = d+ -\-dxl\i = t + + dxjii - d, (6.b.22) 
since t = 2d when s = 1. Again, Tx ~ t when 1. Furthermore, 

(6.b.23) 



d + ^^d^ + xd/ ji 



and one computes 



d^d^ + dx/fi + + g 
" {d + s/^PTd^Y ./d?Tl^ ^ ' 



confirming that Vx is an increasing function of x when s — 1. 

As an example, ior m = g with d = 0.07, K ~ 5.25, s = 1, and with a = 5 = 0.05, rf = 1.5, 
and S — 10/10, Model 2 has an optimal return of 18.7 for the strategy x = 1 with equitable partition of all 
resources amongst 7 males, while Model 1 has an optimal return of 18.4 for the strategy x — 0.3 — l — 7{a + d) 
by investing 30% of resources equitably between 2 males and 70% equitably amongst 7 females. The return 
in Model 1 from the strategy a; = 1 is only 11.4. 

6.C Sigmoid Males 

The theoretical discussion for this case is essentially identical to that for DMR males. Note that the 
consideration of tjjx is really only relevant when x > p (since Ux — I for x e {d,p]) where m" < as for DMR 
return fimctions. Thus, Models 2 & 1 may be solved as indicated in §6.b. 

The linear approximation to Model 2, however, is another matter. Recall from §3.e that for sigmoid 
males M{x) = m{x), on [0,p]. The linear approximation of by A4"'{x) = {x/t)m{t) cannot be expected to 
provide a good approximation to m on [0,p]. The larger is j>, the more pronounced will be the disagreement 
between Ad and Of course, the smaller is p the more a sigmoid approximates to a DMR function. 

For case studies of sigmoid return functions, I employed a sigmoid return function based on the logistic. 
Consider, for x > d, 

ZL 

Kx) ■■= 77 TT^ 7-^> ^ - Z, (6.C.1) 
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with L> Z >0, then l{x) 



L — Z as X ^ oo and l{d) = 0. The point p of inflexion of I is found to satisiy 



r= ^ 6.C.2 
p — ct 

Through this equation, L, Z, and p can be taken to be the model parameters (specifying p provides more 
direct control of whether the sigmoid is more like an IMR or DMR return function, or balanced between 
these two extremes). Substituting (6.C.2) into (6.C.1) yields 

ZL 

= Z. (6.C.3) 



iL-z)[T:^]-+z 



One may compute that 



for the maximum rate of change of I, and 



For I to display its asymptotic behaviour on the interval {d, 1] , one requires r large, which can be achieved 
by choosing Z small relative to p ~ d. 

Substituting (6.C.1) into (3.4.7) yields 

rite'^^*-'*) = Ze^''^*-'') + {L - 2Z)e''^*-^'^ - (L - Z), 

which is not analytically solvable for t, though one can numerically approximate t for given parameters by 
plotting m' {x)x — m{x) as a function of x to locate the zero. Thus, the optimal n for any x will have to be 
determined numerically, either by computing nm{x/n) for the allowed values (A. 2) of n or by numerically 
approximating t and using (3.d.8). 

For this return function, the linear approximation to A4 was indeed poor on [0,p]. Examination of the 
function for (6.C.1) indicates it is much more peaked than for g and h (and does become concave up 
towards the right-hand end of its domain) so the approximation of taking Ux to be v^, = x/t is less accurate. 

The difference between sigmoid and DMR males is most pronounc;c;d, of course, when p is a significant 
fraction of the available resources. With L = 1Q , Z = 0.01, d = 0.05, p = 0.45, a = 5 = 0.05, r/ = 1.5 then 
w = 0.1. With 5 = 1, Model 2 has an optimal return at x = 0.6 = 1 — 4(q! + 5) (this value indicates the poor 
performance of the linear approximation) achieved by investing 60% of resources in a single male and the 
remaining 40% equitably amongst 4 females. Model 1 achieves its optimal return with the same strategy, 
independently of the values of (f) and ji (subject to (j)/ ji = 1). If 5 = 4, the optimal strategy becomes x = 1 
for Model 2 and for Model 1 with large populations. When ^ = 6 and (j) = 24, the optimal strategy for Model 
1 switches to a: = 0.6 again (investing in one male as before), and remains so for smaller values of ^. This 
relative stability of the optimal strategy in Model 1 relative to the value of /z in a fixed ratio </)//x reflects 
the fact that one must invest more than p in each male to obtain a useful return and here p is a substantial 
proportion of the available resources, thus severely limiting the number of males one can invest in. 



7. STEP-FUNCTION MALES AND NON-STEP-FUNCTION 
FEMALES 

The male return function m is given by (3.a.l), whence M. is given by (3. a. 2). The packaging cost for 
females is 5; assume that a + d + 5 < 1, equivalently a + < 1 — ^, so that investment in both genders is 
possible. 
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Results for Model 2 can be obtained from the results for Model 2 in §5.2 as follows. One has 



R^{x) = M{x)S + T{\ -x) = n^TmS + T{1 - x), 

where = [x/{a + d)],. Substitute z = 1 — x and multiply by 1/S to obtain 

R2{z)=J^iz)S + k,rm, (7.1) 

where S := 1/S and kz = rix = [(1 — z)/{a + d)]i. By interchanging the roles of males and females, the 
optimal strategy for (5.3.1) can be obtained from the appropriate case in §5.2. In particular, with 

1 = ^{a + d) + e V & nonnegative integer < e < a + d, (7.2) 

and choose t so that oj := £{a + d) < 1 — S but {£+l){a + d) > 1 — 6, analogous to (6.1-2). Then the optimal 
return for R2, whence for R2, must occur at one of 

0, a + d, 2{a + d),...,£{a + d) = w, 1. (7.3) 

For example, by this procedure one obtains the analogue of (6. a. 6) for Model 2 with step-function males and 
linear females 

f{z) = (5{z - 5); (7.4) 
namely, the optimal strategy occurs at the x value X = 0, 1 — w, or 1 according as: 

i?2(0) > -R2(l) & /?(a + d) > Sr^, then X = 0; 
^2(1) > -R2(0) & /?(a + d) > 5r„, then X = 1; 
i?2(0) > i?2(l) & /3(a + d) < Sr^, then X = l-uj; (5.3.5) 

(5 < e, then X = 1 — w; 



Ri{l) > i?2(0) & P{a + d) < Svm, & 



S>e, then X = 1, 



where e is as in (7.2). 

Due to the asymmetry between males and females in Model 1, however, interchanging the roles of males 
and females does not convert the current cases into those of §6. Since, for the return functions of §3, optimal 
investment always requires equitbale partition of resources amongst several individuals of the same gender, 
i?i can be obtained by maximizing 

——{4> + k) + kf — — 
ji + n \ 

with respect to n and k. As n/ {fi + n) is an increasing function of n, one gets 
where = [x/{a + d)], and 

R,{x) = ^^{cp + fci_.) + (\^) , (7.7) 

where ki-x optimizes (7.6). Now, Cx ■= nxrm/{^J^ + rix) is constant on subintervals of the form 

[0, a + d), [{a + d),2{a + d)),..., [i{a + d),l-6) (7.8) 
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where £ is as above. So, if x < y within one of these subintervals, then 



by the optimality of ki-x 
since / is an increasing function 



i.e., i?i is nonincrcasing on each subintcrval of (7.8), while on [1 — S, 1] i?i takes its maximal value at 1. 
Hence, the optimal return is achieved at one of x values of (7.3) after all. 

Thus, for step-function males, the optimal return for both Models 1 & 2 must occur at one of the x 
values listed in (7.3) and to this extent the situation for both models is analogous to §6. While the analogy 
is complete for Model 2, to underscore the difference for Model 1 between §§6 and 7, recall that for linear 
males and step-function females, the optimal strategy involved investment in at most one male. Consider 
now Model 1 for step-function males and linear females. With notation as above, on [1 — ^, 1], Ri{x) = Cx4> 
while for z > 5 (i.e., x < 1 — 5) one must optimize 

Cx(</' + k) + (3{z - k5) = Cx(t) + (Cx - PS)k + (3z, 

with respect to k, subject to fc < z/S. On [0, a + d), Cx = and the coefficient of k is negative. Whenever 
the coefficient of k is negative, the optimal value of is 1. But when the coefficient is positive, the optimal 
value of k is the largest possible value subject to the constraint k < z/6. The coefficient of k is positive 
iff Cx > i.e., iff the net return from the investment in males per female in the post-release population 
exceeds the total cost due to female packaging cost per translocated female. If so, then it pays to invest in 
more than one linear female, whence the case of step-function males and linear females differs in this respect 
from that of linear males and step-function females. 

Optimization of (7.6) with DMR or Sigmoid females to obtain Ri will be considered in the more general 
setting of §9. 



Ri{x) = Cx{(t> + ki-x) + ki-xf 

> Cx{^ + kl-y) + kl-yf 
>Cy{4>+ kl-y) + kl-yf 

= Riiy), 



l-x 



ki- 


-X 




X 


ki- 


-y 


1 - 


y 


ki- 


■y . 



8. LINEAR MALES AND FEMALES 

Suppose that the individual return ftmctions for males and females are both linear: 

m(x) = l^^ ifa;e[0,d]; f(.) = |0' if [0,5] 

\aix-d), ffxe[d,l]; ^^^> \aiz-5), if [,5,1]. ' 

From (3.b.3), ^A{x) = m(x) and J^{z) = f{z). li d + 5 > 1, then one can only invest in either males or 
females. The returns for Model 1 will be a(l — d)(j)/{ii + 1) and a{l — 5) respectively. 
Now suppose d-\- 5 < 1, i.e., d<l — 5. Consider first Model 2. Then (4.2) is 

{a{l- X- 5), \ix £ [0,(i]; 
V'(x), a X e[d,l ~ 5]; (8.2) 

a{x-d)S, if X & [l - 5,1]; 

where 

V' := a{x - d)S + a{l - x - 6) = {aS - a)x + a{l -6)- adS. (8.3) 

Note that i?2 is continuous and piecewise linear. i?2 decreases linearly on [0, d], increase linearly on [1 — (5, 1], 
and either increases or decreases linearly (as tp{x)) on [d, 1 — 6]. Thus, there are no local maxima on (0, 1) 
and the maximum of i?2 occurs at x = or 1. Hence, if 

QjS 1 S 

i?2(l) = a{l - d)S > a{l -5) = R2{0), equivalently — > (8.4a) 
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choose X = 1, i.e., assign all resources to a single male; while if 



i?2(l) = a(l - d)S < a(l -6)= i?2(0), equivalently ^ < (8.46) 

choose a; = 0, i.e., assign all resources to a single female. Thus, a very simple strategy results. 

Now consider Model 1. On [0,d], M{x;n;xi, . . . ,Xn) = and when investing purely in females the 
optimal strategy is to invest in a single individual, so Ri{x) = R2(x) = Rsix) on [0, d]. On [1 — (5,1], 
F{z; k; Zi, . . . , Zk) =0, so (4.1) reduces to M(x; n; xi, . . . , Xn)Sp = a{x — nd)[<p/ (/i + n)], 1 < n < x/d. This 
expression is maximized for n = 1, so i?i(x) ~ a{x — d)(j)/{fjb + 1) = [m/(m + ^)]R2{x) = R^ix) on [1 — 5, 1]. 

On (d, 1 — 6), one must maximize the expression 

a(x-nd)- \-a(l-x-k5) < n < -, OA; < — —, 

noting that both n = and A: = are non optimal. The first summand is maximized with respect to n by 
choosing n = 1, yielding 

a{x-d)'^^+a{l-x-k8), l<k<^--^, (8.5) 
M + 1 



which is linear in k. The coefficient of k is 

a{x — d.) 



- a5, (8.6) 



which is an increasing linear function of x on (d, 1 — 5), initially negative and with zero at 

X* = ^6{n + l) + d>d. (8.7) 

When (8.6) is negative, (8.5) is maximized by choosing k — 1. But when (8.6) is positive, (8.5) is maximized 
by choosing the largest value of k subject to fc < (1 — x)/S, , i.e., the largest k satisfying x < 1 — k6; so if 
X G [1 - (^ + 1)6, l-e6), choose k = e. 

For simplicity, first suppose that x* > 1 — 6, i.e., 

a l-d-6 



in which case (8.5) is maximized by choosing /c = 1 at each x E {d, 1 — 6) and Model 1 reduces to Model 3. 
Substituting (</)+ l)/(/x+ 1) for 5* in (8.3), the result, denoted ip^, is still linear in x. We have already noted 
that i?3 = i?2 on [0,d] and that R3 = + 1)]-R2(a;) on [1 — i5, 1]. Hence, Ry, is, like R2, piecewise linear 

on [0, 1] and continuous ai x = d. But unlike i?2, R3, is discontinuous at x = 1 — (5, due to the discontinuous 
behaviour of iSp there. When x <l — 5, one invests 1 — x in a single female and this female contributes to Sp. 
The limiting return for iis as x approaches 1 — ^ from the left is lima;^(i_5)- Rz{x) = lima,^(i_5)- ^3(2;) = 
a(l — 5 — d){(f) + l)/(/U + 1) > a(l — 6 — d)cf)/{^ + 1) = i?3(l — 5), i.e., the actual return at 1 — (5 is less 
because one can no longer invest in the female and it drops out of the factor Sp. By the same argument 
as for Model 2, one deduces that there is no local maximum for R^ on (0, 1). Thus, the optimal return for 
i?3 = Ri, granted (8.8), occurs at x = or 1; unless V'3 is increasing, in which case it is possible that the 
limiting strategy lima,_>(i_5)- R^ (x) exceeds both x = and 1. If so, there is a strategy that assigns 1 — 6 — e 
to the male and e + 5 to the female which is more optimal than either x = or x = 1 . Decreasing e increases 
the overall return as long as c > 0. In this case, there would be no actual optimal strategy (i?3 has no 
maximum value on [0, 1]), though one can achieve a return arbitrarily close to the limiting strategy (which 
is the supremum, the least upper bound, of Rz on [0, 1]). This possibility arises only if ^/)3 is increasing, i.e., 
a/a > {fi+ l)/{4> + 1), which is not at all implausible, and if the limiting strategy exceeds the strategies at 
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X = and 1. The limiting strategy will exceed i?i(l) only when 5 < {1 — d)/{l + (j)), which will fail when 
the number of females in the target population is sufficiently large. 

Returning now to the general scenario for Model 1, Ri{x) = Rsix) on [0,x*) with a single female 
invested in on this subintcrval. At x — x* , (8.6) is zero and k drops out of (8.5) and Ri{x*) is independent 
of the value one chooses for k. Suppose x* G [l - {k* + 1)6, 1 - k*d). Specify fc = A:* at x*. 

For X G {x*, 1— (5), (8.6) is positive whence (8.5) is maximized by choosing k as large as possible subject to 
k < {l — x)/S. Denote this value of k by ki^x- For linear returns, one may take the partition of resources 1 — x 
amongst kx-x individuals as equitable partitioning. Observe that on Ij. := \l — {k + 1)5, 1 — k5) fl (a:*, 1 — 5), 
k\-x = k and i?i is linear in x: 

Ri(x) = a(x - rf)^i%i +a(l-x- ki-J). (8.9) 
fi + 1 

Thus, Ri is: piecewise linear on [0,1]; continuous on [0,a:*); but, if k* > 1, discontinuous at x = a;* , 
1 — k*6, 1 — (fc* — 1)6, . . . ,1 — 6. Note that the coefficient of x in (8.9) depends on ki-x and decreases with 
increasing x. If, therefore, this coefficient is negative at x*, then it is negative throughout [x*,l — S) and Ri 
is piecewise linear and decreasing on this subinterval. It also follows that Ri is decreasing linear on {d,x*] 
(since Ri is obtained on this subinterval by substituting 1 for ki^x in (8.9)). In this case, the optimal return 
occurs at a; = or 1. If, on the other hand, the coefficient of x in (8.9) is positive on [x* , 1 — k*6), then Ri 
is linear increasing on this subinterval and either increasing or decreasing on [1 — k*6, 1 — {k* — 1)6), and 
so on. Once Ri becomes decreasing on one of these subintervals, it is decreasing on the remaining such 
subintervals. In this case, the optimal return can only occur at 

a; = 0, a;*, 1 - k*6, l-{k* - 1)6, ...,1-26, 1 (i?i(l) > i?i(l - 6)). (8.10) 

But as we saw in the case when x* > 1 — 6, it may be that Ri does not, because of its discontinuities, have 
a maximum value on [0, 1] , but rather approaches a supremum, achieved as 

lim Ri{x), foiy = l-k*6, 1 - {k* - 1)6, . . . ,1 - 6. (8.11) 

In this case, one must examine the returns at a; = 0, 1 and at any other point of (8.10) for which i?i is 
decreasing on the subinterval that point starts, and the limiting strategy at any point of (8.11) for which 
_Ri is increasing on the subinterval that point terminates. The limit of Ri{x) as x approaches 1 — A;^ is just 
a(l - k6 - d){(j) + k)/{n + 1), which exceeds i?i(l) if and only if 5 < (1 - (i)/((/) + k). 

With d = 0.05, a = 5, 6 = 0.05, a = 1, (/> = 10, /z = 5, one has i?i(l) = 7.92, i?i(0) = 0.95, x* = 0.11, R-l 
is increasing on each of the relevant subintervals Ik and so one must examine the limit of Ri as x approaches 
the righthand endpoint from the left. One finds that as x approaches each of 0.8 and 0.85 from the left, the 
limit of R\{x) is 8.75, this value being the supremum of Ri on [0, 1] (the other limiting values are strictly 
less than 8.75). In this case, there is no optimal strategy, only an optimal limiting strategy: take x close to 
but less than 0.8 (with A: = 4) or close to but less than 0.85 (with k = 3). In practice one can choose an x 
close to either 0.8 or 0.85 that yields a return near 8.75, which is greater than the return of 7.92 at a; = 1. 
Whether such differences in return will prove of practical significance is another question. 

9. DMR/SIGMOID MALES AND FEMALES 

Finally, we suppose that males and females each have either DMR or sigmoid returns (in any combina- 
tion). We may treat these returns together since by the arguments of §3.d-e the optimal strategy for investing 
given resources in a gender with either kind of return is equitbale division of those resources amongst the 
number of individuals characterized as in (3.d.7-8). The packaging costs are d for males and 6 for females 
and we suppose that d + i5 < 1 so that investment in both genders is possible. The male individual return 
function m is either DMR or sigmoid and therefore dtermines a number t satisfying (3.d.7). Similarly, the 
individual female return function / is either DMR or sigmoid and also determines a number satisfying (3.d.7) 
which we denote by T, i.e., T is the unique number satisfying f'{T) = f{T)/T. 
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For Model 2, the overall return function R2 is given by (4.2), which can be written 



R^(x)^n,m(^]s + k,f{^], (9.1) 



where z = 1 — x, Ux is the optimal number of males in which to invest resources x equitably, kz is the 
optimal number of females in which to invest resources z = 1 — x equitably, and it is understood that the 
first summand vanishes on [0,d] and the second on [1 — 6,1]. We resorted to numerical solution of this 
model. For any x, one can find Hx using (3.d.8) when one knows t and kz using the analogue of (3.d.8) for T. 
Alternatively, one can compute nm{x/n) for each n satisfying (A. 2) and select the optimal n and kf{z/k) 
for each k satisfying k6 < z and select the optimal k. Repeating this process for a suitable range of x values 
provides a graph of R2 from which one can determine the x value (to any desired degree of accuracy) yielding 
the optimal strategy. 

But, one can make the following approximations: 

nx^Vx = x/t kz « z/T, (9.2) 

Substituting into (9.1) yields the following approximation to R2: 

({l-x)r{T), ifxe[0,d]; 
R^{x) = I {m'{t)S - f'{T))x + f'{T) =: A{x), if x e (d, 1 - 6); (9.3) 
[xm'{t)S, if XG [l-,5, 1]. 

Note that extending A{x) to all of [0, 1] gives the line joining the two points (0,i?|(0)) and (l,i?f(l)). 
R2 is discontinuous at x = d and 1 — 6 and pieccwisc linear on [0, 1]: on [0,d] it decreases linearly from 
ij«(0) = A{0) to (1 - d)f'{T) < A{d) at x = d; on {d, 1 - 6), R2{x) = A{x) either increases or decreases 
linearly towards ^(1 — 6) > i?2(l ~ <^); on [1 — 6, 1], R2{x) increases linearly to R2{1) = A{1). It follows 
that i?2 achieves its maximal value at either a; = 1 or a; = according as A{x) is increasing or decreasing, 
respectively, i.e., R2 attains its maximum at 

n, iim'{t)s>f{Ty, 

\0, ifm'(i)S</'(T). ' 

Thus, if the approximation (9.3) is effective, Model 2 has the simple solution (9.4), which is to invest all 
resources equitably in males if m'{t)S is greater than f'{T), and equitably in females otherwise. Recall from 
§6.b that the approximation (9.2) never underestimates the returns, so R2{x) < R2{x). 

Numerical examples using the fimction g of (S.d.lO) suggest that the rule (9.4) is not without merit for 
predicting the optimal strategy for Model 2. Further study of the approximation (9.2) may be worthwhile 
for particular applications. That said, (9.4) is not perfect. With S = 1.5, d = 0.2 and m(x) = 6 (1 — [d/x)), 
and 6 = 0.07 and f{z) = (1.5) (l — (t^/z)''), i.e., both m and / of the form g in (S.d.lO), (9.4) predicts 
that the optimal strategy occurs at a; = 1. i?2(l) = 10.8, but, from the graph of R2, the maximal value is 
i?2(0.89) = 11.02, indicating that a slightly higher return is achieved by investing some (0.11) resources in a 
single female. In this example, (9.4) fails because (9.2) overestimates the return at x = 1 {R2{1) = 11.25) 
excessively. Whether such differences in return prove of practical significance is another issue. 

Now consider Model 1. Because resources must be equitably divided amongst individuals of a given sex, 
see §3.d, Ri can be obtained by optimizing the expression 

/x\ f 4) + k\ , , f 1 — x\ „ X „ , 1 - .T . 

where it is understood that n = only on [0, d] where the first summand is understood to vanish and A; = 
only on [1 — i5, 1] where it is understood the second summand vanishes. 

The factor if :=(/> + fc in the first summand of (9.5) merely scales the magnitude of this summand 
and does not affect the optimization of the first summand with respect to n. The first summand in (9.5) is 
therefore of the same form as in (6.b.6). Replacing (6.b.8) with 



il}x{v) := vm 



K 

11 + V 



25 



the analysis proceeds exactly as in §6.b with the result that has a unique relative maximum on its domain 
(0,a;/(i), which is its actual maximum, occuring at Vx = x/tx, where Tx solves (6.b.l2). Hence, the optimal 
value Ux of n is the one of the two integers closest to yielding the greatest value of the first summand in 
(5.5.2). Of course, it is always the case that n^; = on [0, d\ and = 1 on {d, 2c?]. After optimization with 
respect to n, (9.5) becomes 

bx{4> + k) + kf(^^\, (9.6) 



where 



k 



hx:=m[^][-^], (9.7) 



H + Ux 



It remains to optimize (9.6) with respect to k. On [0,d], is zero and (9.6) reduces to fc/((l — x)/k), 
which is optimized as in (3.d.7-8), with / and T taking the places of ni and t respectively. For x € [1 — ^, 1], 
A; = and Ri{x) = hx4>- On (d, 1 — ^), to determine Ri it suffices to optimize 

hxk + kf(^^^, (9.8) 

noting hx > Define 

ix{w) := bxW + wf (^—^) , (9.9) 



w 



on {0,{l-x)/6). Then, 



a^) = b. + f(^—^)-r(^—^) (^—^) (9.10a) 



w / \ w / \ w 



CH = r{'—^)(^-^-^]. (9.10b) 



w 

Thus, £^x is concave down on its domain, and since ^x is assumed smooth, if there is a relative maximum for 
^x it is unique and also the maximum of £,x on its domain. Suppose then that ^x{wx) = and define 

cTx := ^ (9.11) 

Wx 

so that, from (9.10a), 

bx = f'((ya:)^x - f((y.) (9.12) 

is the defining equation for Gx- Since hx > 0, one requires f'{cFx) > f{o'x)/o'xy whence, from the shape of the 
graph of /, ax < T, where as above T is the solution of (3.d.7) for /, i.e., 

S <ax<T. (9.13) 

Thus, Wx ~ {1 — x)/(Tx > (1 — x)/T, i.e., in Model 1 one expects to invest in more females than one would 
just to optimize the second summand of Ri. 

If / is of the form h = k{x - Sy in (S.d.lO), (9.12) gives 

bx = Kiax-5y-'{ax{s-l) + 6), (9.14) 

which is not explicitly solvable for ax- On the other hand, if / is of the form g = k{1 — (S/xY) in (3.d.l0), 
(9.12) yields 

ax = ^/ITlSfl-^ = T^l-^ < T. (9.15) 

y bx + K y bx + K 
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Thus, for functions of the form g, one obtains an expHcit expression for a^- To determine the vahie of k 
maximizing (9.8), one need only test the two integers closest to ax and choose the one yielding the largest 
value of (9.8). 

If therefore, m is of the form h, ot g with s = 1, and / is of the form g, one can compute t, T, Tx and cr^ 
explicitly, which reduces the amount of numerical analysis required to determine Ri. Whatever the forms 
of m and /, however, for any value of x only finitely many numerical computations are required to compute 
i?i(x); specifically, for a given value of x, compute the first summand of (9.5) for the allowed values of n 
and thereby determine Ux, then compute (9.5), with Ux substituted for n, for the allowed values of k and 
thereby determine kx, the optimal value of k, whence i?i(x). 

Numerical studies using the DMR functions of (3.d.l0) with various model parameters revealed no 
obvious restrictions on the value of x yielding the optimal return. To give but one example, using the 
function g of (S.d.lO) for both males and females, with m{x) — 6(1 — (0.1/x)) and f{z) = (1.5)(1— (0.1/z)'^), 
and with /i = 5, the optimal strategy occurs: at x = 1 when </> > 19, with investment in 4 males; at 
approximately x = 0.87 when 16 < </> < 18, with investment in 3 males and one female; at values of x in 
the interval (0.73, 0.77) when 9 < (/> < 15, with investment in 3 males and two females; at values of x in the 
interval (0.44, 0.50) when 3 < < 8, with investment in two males and 4 females; and at values of x in the 
interval (0.24,0.28) when < (p < 2, with investment in one male and 5 females. 

I present one example which illustrates features typified by sigmoid returns. Take both male and female 
returns to be sigmoid. Take Z = 0.01 in (6.C.3) for both return functions, but L = 10 and p = 0.3 for m and 
L = 2 and p = 0.25 for /. Also set d = 0.05 and 5 = 0.1 (packaging costs for males and females respectively). 
With 5=1, model 2 has its optimal strategy at x = 1 with investment in 3 males. As x increases from 
zero to one, Hx increases monotonically from zero to 3, while kx decreases monotonically from 3 to zero. In 
Model 1 (with = /i = 10 and the same model parameter values), there is no optimal strategy but rather a 
limiting optimal strategy as x approaches 0.8 from the left of investing in 2 males and 2 females. Moreover, 
as X increases from zero to one, nx increases monotonically from zero to two but kx no longer decreases 
monotonically. Once investment in a male exceeds its p value, the return on that male apparently affects 
the optimal number of females through Sp in this interesting fashion. 

11. DISCUSSION 

In §§5-10, I have explored at some length various instances of the simple model I have proposed in §§2- 
4. This mathematical study provides solutions of these scenarios, sometimes in the form of a mathematical 
result and sometimes just as a numerical procedure, thereby demonstrating the model is tractable. 

The biological content of the model is primarily to be found in the representation of reproductive return 
of an individual as a function of investment, for each gender, with the further assumption that mS and /, 
where 5 is a suitable measure of mating opportunities of males in the target population, are comparable 
measures. This assumption is taken over from the literature on sex allocation and parental investment 
theory, e.g., Frank (1987). I have assumed that the return functions are nondecreasing functions, initially 
zero to allow for packaging costs, and, with the exception of step-function returns, continuous. Under these 
assumptions, it was found that optimal investment of resources in a given gender was never found to involve 
unequal division of resources amongst several individuals, so that the optimal strategy always amounted to 
equitable partition of resources (though for linear and IMR return functions, the optimal strategy is to invest 
in a single individual). Linklater's (2003) specific proposal of basing a model on the TW effect could be 
achieved by the appropriate choice of return functions, perhaps sigmoid returns for males and DMR returns 
for females in the context of unlimited resources and detailed knowledge of the return on investment of the 
species in question. Limited resources and incomplete knowledge of returns may entail that other return 
functions considered in this paper prove practical for modelling purposes. 

The return functions studied here are mathematical idealizations of likely measures of reproductive 
performance of real organisms. For mammals, say, the return functions might represent the number of 
offspring of an individual of a specific gender in an idealized test population which reach sexual maturity. 
How well continuous return functions serve to represent discrete returns might depend on the biology of the 
species in question and the response to levels of investment. Perhaps multi-step- function returns will be 
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more appropriate. One has some indieation of the outcomes of the model with muhi-step-function returns 
to the extent that they approximate hnear, DMR, or sigmoid returns. In general, however, one must resort 
to numerical analysis to solve the model with multi-step-function returns, beginning with the essentially 
combinatorial problem of determining the optimal division of given resources amongst a given number of 
individuals of the same gender. The solution to this problem will depend on the exact nature of the multi-step 
function and the optimal division of resources need not be equal division between individuals. 

Practical applications require knowledge of the return functions, presumably by experimentation. For 
endangered species, such would require possibly unjustifiable risks. For black rhino, translocations have not 
been carried out with varying degrees of investment under experimental conditions. Rather, translocations 
are typically carried out under a regime of fixed investment. It may be appropriate to model such translo- 
cations by step-function males and females. Survival and fecundity estimates of post-release animals could 
be converted into, say, mean number of offspring/per year (where for males this mean must include inverse 
weighting by its mating opportunities), and this measure employed as Vm for males and rj- for females in the 
modelling performed in §5 for application to future translocations under the same investment regime. Even 
such a rudimentary model may be difficult to parametrize however; male fecundity for back rhino is unlikely 
to be well known since paternity of offspring is rarely known. 

Note that the future reproductive success of an individual is dependent upon its current age. I do not 
model this dependency but implicitly assume that the individuals available for translocation are of optimal 
age for future reproductive returns. In practice, it may or may not be feasible to meet this assumption. If not 
met in a known way, it may or may not be straightforward to modify the model to allow for the particular 
situation. 

The fact that in some circumstances the total return function Ri of Model 1 does not a have a maximal 
value on the interval [0, 1] but only a supremum (least upper bound) is disconcerting. The phenomenon 
arises from the discontinuity in Ri stemming from the numerator of Sp. In particular, an arbitrarily small 
investment in a female beyond the packaging costs makes a contribution to i?i through its contribution to 
Sp. When this investment in that female decreases to an investment at the level of packaging costs, that 
investment will optimally be redirected to another female and the decrease in the number of females invested 
in causes Sp to change discontinuously. As first noted in §8, it may then happen that the value of i?i at 
the point of discontinuity y is less than the one-sided limit of Ri as x increases towards the discontinuity 
y, whence Ri{y — e) > Ri{y) for arbitrarily small positive e. The limiting return lim^.^^- Ri{x) may turn 
out to be the supremum of -Ri on [0, 1], but does not represent a strategy which can be implemented. The 
best one can do is to choose the optimal strategy for x = y — e, for some suitable and small e. It remains 
to be seen whether this phenomena will prove of practical import. Theoretically, however, it perhaps raises 
the issue of whether it is realistic to suppose that returns on female investment increase continuously from 
zero when investment exceeds the known packaging costs. Suppose instead, as with the (multi)-step-fimction 
return, a further threshold level of investment beyond the known packaging costs is required for a nonzero 
return. If at a; = y some female receives exactly the threshold level of investment, as x increases toward 
y, Ri is expected to vary continuously and therefore now take its value at x = y given by this limit, i.e., 
Ri{y) = liui^^y- Ri[x). When x exceeds y, the female in question no longer yields a return and the 
optimal strategy for a; > y is to direct investment to another female already invested in (if such occurs) and 
the numerator of Sp decreases, discontinuously, by one at x ~ y. Now limj,^j,+ ^ Ri{y), reflecting this 
discontinuity at x = y, and one expects this limit to be less than Ri{y). Thus, such threshold values would 
not eliminate discontinuities from R\ but render them innocuous, as in the cases involving step-function 
returns. As with multi-step-function returns, for a return function with a threshold value beyond packaging 
costs for a nonzero return, the argument of Appendix A which prescribes equitable partition of resources in 
a given gender as optimal will be violated. 

Recall that the model of this paper was designed to optimize the returns of the translocated animals 
in line with Linklater's proposal to optimize genetic flow into the target population. As such, the optimal 
strategy provided by this model is not necessarily the optimal strategy as regards reproductive performance 
of the target population, as this latter goal may not coincide with that of maximizing the reproductive 
performance of the translocated individuals. As indicated in the introduction, meta-population management 
involves a number of issues, of which I have considered just one. To indicate one aspect of this point, consider 
a scenario in which male returns are sigmoid, with Z = 0.01, L = 10, p = 0.3, and d = 0.05 in (6.C.3), and 
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female returns are the DMR function /(z) = (1.5)[1 - [O-l/zf] of (S.d.lO). For Model 2, with S = 1, the 
optimal strategy is x = 1 with investment in three males. For Model 1 {S = 1) with /j, > 14, the optimal 
strategy is still x = 1 but with investment in 8 males. As jj, decreases, so does the value of x at which the 
optimal strategy occurs. When /i = 1, the optimal strategy occurs at a; = 0.27 with investment in two males 
and 7 females. These results are neither surprising nor objectionable. When the number of females in the 
target population is small, the predictions of the model may serve the population reproductive performance 
less well. Consider the extreme case when (j) ~ 0, whence S = (other model parameters as above). For 
Model 2, the optimal strategy is x = with investment in 6 females. For Model 1, however, the optimal 
strategy depends upon the value of /x. For /z > 62, Model 1 yields the same result as Model 2, but when 
/i = 61 the optimal strategy occurs at a; = 0.1, with investment in one male and 6 females; some investment 
has been diverted from the 6 females to a single male because, according to Model 1, that investment will, 
through a male's greater reproductive potential, increase the overall return from the translocated individuals, 
despite the presence of the 61 males in the target population. Perhaps more surprising, there are large values 
of fx, beginning with fj. = 44, for which, according to Model 1, the optimal strategy favours investment in 
males at the expense of the number of females translocated, e.g., for /x = 20 the optimal strategy occurs at 
X — 0.36 with investment in 3 males and five females for a net return of 10.1 versus a return of 7.1 for the 
strategy of a; = 0. Unless there are extenuating circumstances regarding the genetic constitution of the males 
of the target population, in a case when there are few females in the target population and Model 1 predicts 
an optimal strategy favoring investment in males at the expense of the number of females translocated, other 
considerations may veto the recommendation of this strategy. As remarked previously, however, it remains 
to be seen whether such circumstances will prove of practical importance, e.g., differences in return between 
certain strategies may prove of negligible importance in practice. 
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In this appendix, I present a mathematical result useful for obtaining the optimal division of resources 
X invested in a single gender. I use the notation appropriate for males. 

Let m denote an individual return function, with packaging costs d. For fixed resources x, let Vnix) 
denote the set of partitions of x into n amounts. A particular partition ^» = x, Xi G [0, 1], may be 

denoted P{x;xi, . . . ,a;„) e Vn{x). Write yj := xj — d, i = 1, . . . ,n. Since it is wasteful to invest no more 
than d in an individual, one may suppose that Xj > d, i.e., yj > 0, and 
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(A.1) 



whence 



nd < X 



(A.2) 
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which constrains the possible values of n. Thus, yj G (0, x — nd). Put 

n 

M{x; n; yi, . . . , y„) := ^ m{yj + d). (A.3) 
i=i 

One can perform the optimization of M{x; n;xi,. . . , x„) with respect to . . . , a;„ and n as follows. 

First, fix n subject to the constraint (A. 2). The problem is then to first choose the yj so as to maximize 
M(x:n;yi, . . . ,yn), and then to optimize over the allowed values of n. Assuming M{x;n:yi, . . . ,yn), without 
constraints, is a function of (t/i, . . . ,yn), each yj G {0,x — nd), then we have a standard optimization 
problem with constraints. Put h{yi, . . . , = J2]=i 2/j -{x- nd) and L{yi, ...,y„):= M{x; n;yi, . . . ,yn) - 
Xh{yi, . . . ,yn), where A is a Lagrange multiplier. Then, 

— = m'iyj +d)-X j = l,...,n 'dX ^ ^ 

Setting these equations to zero gives the solutions for the constrained optimization: A = m'{yj + d), j = 
1, . . . , n and h = 0, i.e., 

n 

m'{yj + d) = m'{yk + d) for all j, k yj = x — nd. (-^-S) 

The Hessian H oi L is the diagonal matrix 

H = diag(m"(t/i +d),...,m"{yn+ d)) , (A.6) 

evaluated at the solution of (A.5), which is nondegenerate provided each of these second derivatives is 
nonzero. A sufficient condition to guarantee that the solution is a relative maximum is that m" < at the 
solution of (A.5). 

As a function of (t/i, . . . , j/„), the domain of interest of M{x; n;yi,. . . , yn) is % € (0, x — nd), j = 1, . . . , n, 
with the constraint h = restricting M to the hyperplane Vj = x — nd. If D denotes the domain of the 

constrained function, i.e., D is the hyperplane slice of the open domain of M, then its closure is compact and 
at any point of its boundary, D\D, some yj must be zero. As a constrained function of (yi, . . . , y„), for fixed 
X and n, if M has a maximum on its domain D, it must occur at a critical point, which are characterized 
by (A.5). Otherwise, the 'maximum ' of M occurs on the boundary of D. Now, the constrained function is 
continuous, so this maximum can be obtained in a continuous fashion. But, at that boundary point, some 
yj = 0, i.e., one is assigning resources d to the j'th individual. This tactic is nonoptimal. It is no worse to 
reassign the resources d to the other n — 1 individuals, i.e., there is a partition P(x; xi, . . . , Xn-i) € Vn-i{x) 
that is at least as optimal as any partition in 'Pn{x). Note that the return fimction Af(x; n; xi, . . . , x„), now 
viewed as a function of (n, xi, . . . ,x„), may change discontinuously with the change from n to n — 1 in the 
reassignment of the resources d from the abandoned individual. 

The observations of the previous paragraph entail that the optimal investment of resources x occurs at 
a solution of (A.5) for one of the allowed values of n given by (A. 2). In practice, one can begin with the 
largest allowed value of n, solve (A.5) for local maxima, and then proceed to smaller values of n, stopping 
with n = 1 (investing in no individuals is nonoptimal!). The optimal strategy can be identified from this list 
of possibilities (if there are no critical points, the optimal strategy is to invest in a single individual). 



APPENDIX B 

We have seen that for the continuous return fmictions of §3 the optimal investment of resources x 
in a specific gender (say males) may always be construed as equitable division of resources amongst some 
number Ux of individuals satisfying (A. 2): M{x) = nxm{x/nx). Now M{x) = m(x) when x € [0, 2(i], with 
= on [0, d\ and = 1 on (d, 2d\. On (srf, (s + l)rf] , it is possible to invest in s males, and the optimal 
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strategy at x E (srf, (,s + is determined by the relative values of m.{x), 2m(x/2), . . . , sm{x/s). Whichever 
of these functions is largest at such an x is also, by continuity, largest near that x. Indeed, one may view 
the construction of M{x) as follows. Construct the graphs of each of the functions := nm{x/n), for 

those positive integers n satisfying nd < 1. This is a finite number of continuous graphs. On [0, 2(i] all but 
fi{x) are zero. So M{x) will be continued by fi{x) = m{x) beyond x = 2d until one of fn{x) equals and 
then exceeds m{x). This process repeats for the remainder of the domain [0, 1] of M. whence M is seen to 
be made up of pieces of the graphs of fn{x), joined continuously together. Thus, A4{x) is indeed continuous. 

For linear returns m(x) = a{x — d), m„{x) = a{x — nd) which is parallel to m. Thus, Ai{x) = m{x) 
as deduced in §3b. For an IMR return m, observe that m'^{x) = m'{x/n) < m'{x) by definition of IMR 
returns. Thus, as with linear returns, no m„ can intersect m and again M{x) = to(x) as deduced in §3c. 
This same argument, applied to sigmoid returns, establishes that no m„ can intersect m for x < p, whence 
M.{x) = m{x) on [0,p] as argued in §3e. 

DMR return functions give a particularly simple picture. With j > i > 1 and jd < 1, suppose that 
fj{x) = fi{x), for some point Xij G {jd, 1]. Since fj{x) = m'{x/j) > m'{x/i) = fl{x), then at the intersection 
Xij, fj does indeed cross over the graph of fi. If, for fixed i, Xij < Xi^k whenever k > j, then 



M{x) = < 



' m{x), X G [0,a;i2]; 
f2{x), a; e [a;i,2,a;2,3]; 
fsix), a; e [a;2,3,a;3,4]; 



For the functions h and g of (3.d.l0), and with r = 1/s > 1 for h, one finds, respectively. 



Xi.'i — 



-i'')d 



if 



d = ^f + f-'^i H h jis-i + d. (B.l) 



Note that if — ji^ > for r > 1, and both the functions in (B.l) are increasing functions of j for fixed i 
when j > i. Curiously, both functions are proportional to d, but this fact is a coincidence (one can construct 
other DMR functions that do not share this property). Putting s = 1 in the formula for g, one obtains the 
particularly simple result that Xij = {j + i)d. Since 2:^+1,1+2 — Xij^i^i = 2d, in this case 



M = < 



' m(x), a; G [0,3d]; 
2m(a;/2), x€[3d,5d]; 
3m{x/3), xG[5d,7d]; 



(B.2) 



For the sigmoid function of (6.C.1), one does not obtain an explicit formula for Xij. 
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